Evaluating holonomic quantum computation: beyond adiabatic limitation 
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The proposal of the optical scheme for holonomic quantum 
computation is evaluated based on dynamical resolution to 
the system beyond adiabatic limitation. The time-dependent 
Schrodinger equation is exactly solved by virtue of the crank- 
ing representation and gauge transformation approach. Be- 
sides providing rigorous confirmation to holonomies of the 
geometrical prediction that holds for the ideally adiabatic 
situation, the dynamical resolution enables one to evaluate 
elaborately the amplitude of the nonadiabatic deviation, so 
that the errors induced to the quantum computation can be 
explicitly estimated. 



The recently proposed holonomic approach to quan- 
tum computation |Q- |^ surely predicts a striking contri- 
bution to the application of quantum physics. Transcend- 
ing the traditional dynamical means for quantum compu- 
tation, the holonomic approach realizes quantum infor- 
mation processing by endowing the quantum code with a 
non-trivial global topology (a gauge field potential) and 
the associated holonomies then allow for the universal 
quantum computing. Specifically, in the scheme of holo- 
nomic quantum computation, information is encoded in a 
degenerate eigenspace of the governing Hamiltonian and 
the holonomies (abelian as well as non- abehan) |]- § 
are acquired by driving the system to undergo appropri- 
ate loops in the parameter space adiabatically. Besides 
suggesting an intriguing connection between the gauge 
fields and the information processing, such a geometri- 
cal means for quantum manipulation is believed to have 
built-in fault-tolerant features ||]- due to its inher- 
ent stability against local perturbations. Considerable 
attention has been addressed to this topic recently and 
the all-geometrical implementation for universal quan- 
tum gates has been proposed by optical schemes, based 
on laser manipulation of ions confined in a Paul trap jl^ j 
or neutral atoms in an optical resonator p3| . 

The existent exploration for holonomic quantum com- 
putation is based on the analysis by pure geometrical 
fashion. It is true that in the adiabatic limit the holon- 
omy associated with the evolving loop is determined by 
the path traced by the time-evolution ray and the cur- 
vature of the ray space. This involves the abelian holon- 
omy (the Berry phase) and the non-abelian one merely 
known as adiabatic connection. Nevertheless, as a whole 
physical problem, as the dynamics of a system generates 
a time-dependent physical state, a specified geometrical 
object (the ray) is generated as well. In such a sense. 



dynamics determines the holonomy through determining 
the ray itself and its path. Moreover, in view that the re- 
alistic evolution of a physical system could not be ideally 
adiabatic and the nonadiabaticity shall alter the time- 
evolution of the ray and thus inevitably induce deviation 
from the adiabatic consequences. The evaluation of such 
deviation and the resulting errors in quantum computa- 
tion is definitely a dynamical problem that goes beyond 
the geometrical exploration. 

In this paper we employ a tractable model of the op- 
tical scheme to exploit this subject. For the appropri- 
ately chosen loops of the Hamiltonian in the parameter 
space, the time-dependent Schrodinger equation is ex- 
actly solved by virtue of the cranking representation and 
gauge transformation approach. The derived dynami- 
cal evolution of the system recovers the holonomic trans- 
formation provided by geometrical consequences, includ- 
ing the simple abelian phase factor and the general non- 
abclian operation. Thus our results provide further con- 
firmation of the geometrical prediction, and besides, the 
errors caused by nonadiabatic effects for the holonomic 
quantum gate operation can be estimated explicitly. 

For the proposed optical scheme of holonomy quantum 
computation |jl^,|l3| , the basic idea relies on the adiabatic 
passage via the dark states since the dynamical evolution 
restricted to such a space is completely trivial. The sys- 
tem encoding the qubit is realized by a four-level A-type 
trapped ion (or a similar cavity atom). The three ground 
levels \gi) [i — 1, 2, 3) are highly degenerate and each cou- 
ples to the excited state |e) in a tunable way. The states 
\gi) and \g2) stand for the computational bases |0) and 
respectively, and \g^) is an ancillary level required 
for implementation of gate operations. Such a system 
admits two dark states that have no contribution from 
the excited state. Through changing the Rabi frequen- 
cies and driving the dark states to undergo appropriate 
cyclic evolutions in an adiabatic fashion, the universal 
single-bit gate operations e*"^'^^^^' and e"^'^" can be gen- 
erated due to the global geometry of the bundle of the 
eigenspace of the dark states. 

To evaluate the gate operation e"^'^^^^' from a dynami- 
cal viewpoint, let us explore the state evolution generated 
by the periodic Hamiltonian irUTT 



H{t) = n sin e{a2e + (Te2) + ^ COS e{az^e'f + (Tese"''^), 

(1) 

where is a fixed parameter and ip is assumed to rotate 
at a constant frequency 7 for convenience. The equation 
of motion for the system is 
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i^^\^{t))^Hit)\^it)). 



(2) 



It is known that in the adiabatic hmit, the geometri- 
cal exploration shows that the dark state of the system, 
I !?(<)) = cos 0|g2) — sin0e*^*|(73), shall acquire a net Berry 
phase IJI after a period T = 27r/7: \D{T)) = e"l'\D{0)) 
with (j) = 47rsin^ 9. To give a dynamical resolution to the 
system beyond adiabatic limitation, we use the crank- 
ing representation, that is, the Hamiltonian can be 
regarded as a cranked one 



H{t) 



■i'rto-33 



(3) 



where 



Hq = ilsm0{a2e + 0-e2) +^lcOS0{a3e + Ces), (4) 

and the unitary transformation e^^^"^^ can be viewed as 
an element of the SU(3) group. Consequently, the dy- 
namical invariant of the system can be shown as 

I(t) = e'^^*-- {Ho + 7^33)6-^^*"- - H{t) + -fass, (5) 

which satisfies lfl4ll 



dl{t) _ dl{t) 



dt 



dt 



i[I{t),H{t)]^0. 



(6) 



The second term of (^) accounts for an extra gauge po- 
tential since H{t) depends on time explicitly. Now the 
recurrent basis \ip{t)) of the system [it differs from the ba- 
sic solution of the Schrodinger equation (|^) only by 
a phase factor] can be obtained by solving the instanta- 
neous eigensolutions of I{t). It turns out, one only needs 
to solve the characteristic equation 



- ij/n)x^ -x+ {-f/n) sin^ 0^0 



(7) 



and the eigenvalues of I{t) are given by Ei{^) = 
(i = —1,0,1). It is straightforward to show that the 
recurrent basis |V'o(0) represented by the middle number 
Eq approaches to the dark state \D{t)) in the adiabatic 
limit j/il 0. Now the leakage error induced by the 
nonadiabatic effect can be conveniently estimated by the 
overlap [see Fig. 1(a)] 



v{o,^) = \{Mo)\d{qW^\{Mt)\d{t))\^^ 



(8) 



Besides the leakage, the nonadiabatic evolution shall 
result in deviation to the desired phase factor. It fol- 
lows, instead of the net Berry phase, the cyclic evolution 
here induces a total phase (the so-called Lewis-Riesenfeld 
phase) 

$ = / {i:o{m-^-H{t)\Mt))dt = Eo — . (9) 







dt 
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The detailed depiction of the deviation for the phase fac- 
tor is shown in Fig. 1(b). Noting that in the adiabatic 
limit, the total phase [na| 



$ = lim 2n 



7/0 



An sin 6*, 



(10) 



the geometrical consequence is thus recovered. 

The validity of the above evaluation is based on a 
presumption that the initial state |-D(0)) = cos 6*192) — 
sin 6*153) can be generated from the computational ba- 
sis [(72) and so the inverse process. Explicitly, such pro- 
cesses can be accomplished by the driven Hamiltonian (|l|) 
through changing the parameter 9 adiabatically. Con- 
ventionally, the nonadiabatic effect here shall lead to an 
additional error for the quantum computation. However, 
such an error can be in principle avoided through ap- 
pending a matching interaction to compensate the gauge 
potential term induced to the system. Specifically, one 
can use the following Hamiltonian (setting ip — 0) 

Htot{t) = H{t)+Had{t),Had{t) = i9{t)ia23 - (732). 

(11) 

It follows that the dynamical invariant of the system 
Htot{t) now has a form I it) = H{t), thus the above state 
transformation can be processed exactly. Physically, the 
interaction Had{t) can be realized by a microwave cou- 
pling to the two degenerate levels [(72) and [53), with its 
intensity accurately controlled through a derivative feed- 
back process. 

Now we investigate the gate operation e*"^'^" achieved 
by the holonomic means. The corresponding evolution is 
generated by the Hamiltonian 

H{t) = r2sin6'c0S93(CTie + Cel) 

4- sin 6* sin (^((726 + o'e2) + cos 9{a3e + o'es), (12) 

where the parameter ip = ^t. As is known, the adi- 
abatic cyclic evolution of the Hamiltonian generates a 
non-abelian holonomy due to its degeneracy structure of 
the dark states. It can be easily worked out, from the 
formula of Ref . [Q , that the holonomic transformation 



uc 



i27T COS ODjj 



(13) 



where Dy = i{\D2){Di \ - \Di){D2\), and the two dark 
states, \Di) = cos9\gi) — sin 9\g3) and ID2) = [52), span 
the degenerate space of the starting (ending) Hamilto- 
nian. Note that the Hamiltonian (021) possesses an su(4) 
Lie algebraic structure and dynamical resolution to the 
system is usually very complicated. Surprisingly, as we 
shall show in the following, this system can be exactly 
solved by the gauge transformation approach p6| , p7t , 
and its dynamical evolution analytically manifested thus 
leads to a complete understanding of the adiabatic and 
nonadiabatic properties for the time-dependent Hamilto- 
nian system. 

Similar to the cranking method used above, we intro- 
duce the unitary gauge transformation 



-7t(cri2— 0-21) 



(14) 
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to the equation of motion for the system, from which a 
covariant Schrodinger equation is stemmed 



i^^\^g{t))=Hg\^g{t)), 



(15) 



with the gauged Hamiltonian 



iU: 



» dt 
= Q sin 0{(7ie + cTgi ) 



(16) 



In view that the above Hamiltonian is time independent, 
the basic solutions |5'g(i)) to the covariant equation ( p^ ) 
can be easily obtained and the corresponding eigenvalues 
are as follows 



Ei^2 = 



1 -4(^)2 C0S2( 



1 



where 



1 _4(_L)2cos2( 



COS 9 



1/2 



1/2 



1 + {j/ny 



(17) 



(18) 



The dynamical basis of the system (g^ can be directly 
obtained as |^'„(i)) — Ug\'^g{t)), from which one can see 
that En has the natural implication related to the total 
phase. Now the time evolution operator generated by the 
Hamiltonian ( p^ ) can be given 



Uc{T) 
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n=l 
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n=l 



*„(r))(vi/„(o)| 

3-^^"'?|*„(0))(M'„(0)|. 



(19) 



Considering the asymptotic behavior of the evolution 
in the adiabatic limit, it follows that lim^/f2_^Q ~ 
±cos6, and the phase-equipped dynamical bases 1^*1 (t)) 
and \^2{t)) have the form 

\^i{t)) = ^e-*'^*™^''[(cos6'cos7t + isin7t)|5i) 
+ (cos6'sin7i — i cos7t)|g2) — sin^l^a)], 

|*2(i)) - ^e*''*™"^[(cos6'cos7i-isin7t)|5i) 

+ (cos6'sin7i — icos7t)|52) — sinOlg^)]. (20) 

One can verify that they are the instantaneous eigen- 
states of the Hamiltonian ( |l^ with a two-degeneracy 
eigenvalue 0, and the equipped phases are just the Berry 
phases accordingly. Thus the cyclic evolution restricted 



to the space spanned by these two states is purely geo- 
metrical and can be denoted as 



uiT) 



-■1277 COS 6 



|*i(0))(*i(0)Ke 



z27r cos 9 I 



*2(0))(*2(0)| 

(21) 



with |*i(0)) = (|Di)-i|D2))/V2and |*2(0)) = (|i?i) + 
i\D2)) / \/2. It can be easily recognized that the opera- 
tor (]2l|) is just the non-abelian holonomy (|l3|), thus the 
geometrical nature is verified again. 

The above dynamical resolution to the system is im- 
portant. Besides offering a vivid verification to the re- 
markable formula of non-abelian holonomy 0|, which 
holds for the ideally adiabatic situation, it enables one to 
evaluate elaborately the amplitude of the nonadiabaticity 
deviation and the resulting errors to the holonomic gate 
operation e*'^'^" . In detail, the population transfer from 
the initial state |^(0)) — \g2) is pictured in Fig. 2. The 
leakage out of the computational space can be estimated 
by the projection (see also Fig. 2) 



'?(^'S) = Ei(^^i^(^)i*(o))i' 
1=1 



(22) 



Similar to the former case, to transform the computa- 
tional basis \gi) into the dark state \Di) — cos 6*151) — 
sin^l^s) and to invert the process successfully, one needs 
to use the Hamiltonian ( p^ ) (with 9 tunable and (p = 0) 
along with a matching interaction 

Htotit) = H{t) + Hadit), Had{t) - z^(i)(ai3 - ^31)- 

(23) 

Up to now, we have investigated the single-qubit holo- 
nomic operations of the optical scheme and evaluated 
the nonadiabaticity-causing errors for quantum compu- 
tation. It deserves to point out that, the proposal |1J] of 
geometrical implementation for the controlled two-qubit 
phase shift gate, e*"^'^^^^^^', which is sufficient for the uni- 
versal quantum computation along with the two single- 
qubit gates, can be explored in a similar way. The scheme 
is realized by two-color laser manipulation on the trapped 
ions [Q. Briefly, the transition from the ground states 
132) and Iga) to the excited state |e) is driven by two 
different bi-chromatic laser beams with their amplitudes 
and phases of the Rabi frequencies controllable. The fre- 
quencies of the laser fields are tuned so that the two- 
photon process, exciting pair ions, is resonant and the 
single-photon excitation is off-resonant. Hence the sys- 
tem can be described by an effective Hamiltonian in the 
notation of Ref. ]l^ 

H,ff cx -|r!i|2(e*2^i|ee)(5252| +/i-c.) 

+ |r!2p(e'2^^|ee)(5353| + /i.c.), (24) 

where the relative intensity of the Rabi frequencies 
tan0 = — |r2i|2/|r22p and the phase difference (/?/2 = 



3 



(fii — f2 are tunable. One can see that the bases \gigi), 
I5152), and I (7251) are decoupled from the evolution, and 
the component \g2g2) serving as the code |11) evolved in 
an enclosed space spanned by {\g2g2), I.9353), |ee)}. In- 
troducing the su(3) generators explicitly 



Ae2 

A23 



e'2^Mee)(5252|, A,^ ^ e'^^'\ee){g3g3l 



15252) (5353 1 



Al 



A,, 



(25) 



the Hamiltonian can be rewritten as 



Heff = g sine {A2e + Ae2) + g cose {Asee^'f + Ae3e '"^). 

(26) 

Obviously this Hamiltonian possesses an su(3) algebraic 
structure isomorphic to that of system (Q), thus all the 
discussions therein also hold for the present system 



It should be noted that, the effective Hamiltonian (24), 
respecting a second-order process of the interaction, is 
quite a rough description of the model. Specifically, it 
ignores the same second-order process induced by virtual 
photons excitation in the self-transitions of the states: 
|52(3)) |52(3)) and |e) |e). It can be anticipated 
that such self transitions shall dress the energy levels of 
the ions and lift the degeneracy of the ground states, 
which in turn affects the desired gate operation. Detailed 
exploration of this point shall be presented in a future 
report. 
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Caption of Fig. 1: 

Deviation induced by nonadiabaticity for abelian 
holonomy. (a) The overlap 77 for the parameters 9 g 
(0,7r/2) and -//fl G [0,1]; (b) The total phase $ for 
Oe [0,7r] and -f/fle [0,1]. 

Caption of Fig. 2: 

Deviation induced by nonadiabaticity for non-abelian 
holonomy. The initial state is prepared in |I?2)- The 
two solid curves show the results for the population of 
the target state on |Di) and I-D2), as a function of 1 — 
cos 9, respectively. The dashed curve depicts the total 
population 77 on the computational space. Figures (a), 
(b), (c) and (d) correspond to j/il = 0.01,0.2,0.5 and 
0.8, respectively. 
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